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The dark energy, dark matter and baryon densities in the Universe are observed to be similar, with a factor of 
no more than 20 between the largest and smallest densities. We show that this coincidence can be understood 
via superhorizon domains of randomly varying densities when the baryon density at initial collapse of galaxy- 
forming perturbations is determined by anthropic selection. The baryon and dark matter densities are assumed 
to be dependent on random variables 0^ and 6^ according to p^ m <x 6™ and p/, oc 9^, while the effectively 
constant dark energy density is dependent upon a random variable §q according to pg The ratio of the 

baryon density to the dark energy density at initial collapse, rg, and the baryon-to-dark matter ratio, r, are 
then determined purely statistically, with no dependence on the anthropically-preferred baryon density. We 
compute the probability distribution for rg and r and show that the observed values of rg and r can be naturally 
understood within this framework. In particular, for the case a = 2, p = 1 and n = 4, which can be physically 
realized via a combination of axion dark matter, Affleck-Dine baryogenesis and frozen quintessence with a 
potential, the range of rg and r which corresponds to the observed Universe is a quite natural, with a probability 
which is broadly similar to other ranges of rg and r. 
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I. INTRODUCTION 

A fundamental issue for cosmology is the similarity of the 
observed densities of dark energy, dark matter and baryons. 
There have been several proposals which attempt to explain 
either the ratio of the dark energy density to the total matter 
density or the ratio of the baryon density to the dark matter 
density. Combining these proposals can, in principle, account 
for the observed coincidence of all three densities. 

It has been suggested that the dark energy density could 
be explained by anthropic selection |Q1 0] • This is because 
density perturbation growth is damped in a dark energy dom- 
inated Universe. Therefore if there is a variation of the dark 
energy density between domains of size much larger than the 
presently observed Universe, we might be able to understand 
why the dark energy density is close to the total matter density 
today. 

There have also been many proposals for explaining the 
baryon-to-dark matter ratio. A particularly popular approach 
is based on a conserved charge shared between baryons and 
dark matter, which results in similar number densities of dark 
matter particles and baryons. In most cases this requires that 
the dark matter be asymmetric, with mass in the range 1-10 
GeV. 

If dark matter is not related to the baryon number by a con- 
served charge, then we need to consider a different explana- 
tion for the baryon-to-dark matter ratio. One way is to have 
similar physical processes which independently determine the 
baryon and dark matter densities. For example, right-handed 
sneutrino dark matter can be generated by a process similar to 
Affleck-Dine leptogenesis [3]. This approach allows a wider 
range of dark matter particle mass, since there is no direct link 
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between the number densities of the dark matter particles and 
the baryons. A particular class of such models arises if dark 
matter is a thermal WIMP. In this case the observed baryon 
density must be related to thermal WIMP annihilation fl. 
There have been a number of proposals to relate the baryon 
density to the thermal "WIMP miracle" density Ij48[l, either 
by modifying a pre-existing baryon asymmetry |4l |5[ or by 
generating the baryon asymmetry from thermal WIMP anni- 
hilations or decays 1(1 0]. 

Combining any of these proposals for the baryon-to-dark 
matter ratio with anthropic selection of the cosmological con- 
stant could then explain the coincidence of all the observed 
densities. 

Alternatively, we could consider the possibility that all 
three densities have a domain-dependence. For example, a 
domain dependent dark matter density is natural in the case 
of axion dark matter [9], while a domain-dependent baryon 
asymmetry is possible in Affleck-Dine baryogenesis ifToL [Till . 
In fact, there may be a strong motivation for considering a 
domain-dependent baryon asymmetry. An increase in the av- 
erage baryon density at a fixed temperature will cause a much 
larger increase in the baryon density in galaxies ijgl [HI. As 
a result, the observed baryon asymmetry is likely to be rather 
close to the anthropic upper bound on the baryon asymme- 
try. To understand this additional coincidence may then re- 
quire that the baryon asymmetry is anthropically selected via 
a domain-dependence. The baryon-to-dark matter ratio could 
then be either due to a conserved charge or due to domain- 
dependence of the dark matter density. We will consider the 
latter possibility in the following. 

In 111 211 we presented a general analysis of the baryon-to- 
dark matter ratio in the case where both the baryon and dark 
matter densities have a domain-dependence due to random 
variables which correspond to effectively massless fields dur- 
ing inflation. The key assumption of the model is that the 
baryon density when the perturbations responsible for galax- 
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ies initially become non-linear and begin to collapse is the 
dominant anthropic parameter, such that equal values of the 
baryon density at initial collapse correspond to equal proba- 
bility for the evolution of observers in galaxies. In this case 
the baryon-to-dark matter ratio r is purely statistical and is 
completely independent of the baryon density at initial col- 
lapse lfl2ll . It was shown that observed baryon-to-dark matter 
ratio, r w 1/5, is then a completely natural value in this frame- 
work. 

Our previous discussion did not include dark energy. Here 
we consider whether random domains of all three densities 
can naturally account for the observed coincidence. To do this 
we must compute the probability distribution of the dark en- 
ergy, dark matter and baryon densities for a given dependence 
on the random variables. 

In the following we will analyse the probability distribu- 
tion in a completely general way, without specifying the ori- 
gin of the random variables. The simplest possibility is that 
the random variables are fields which are effectively massless 
during inflation and are randomized by quantum fluctuations 
|[TTl[l2Tl . We will show that domain-dependence of the baryon, 
dark matter and dark energy densities, when combined with 
the principle that the baryon density at initial collapse is the 
primary anthropic selection parameter, can naturally account 
for the observed coincidence of all three densities. 

As an explicit example of a domain-dependent dark en- 
ergy density, we will consider a frozen quintessence model 
111311 with a (j)g potential, where the mass and couplings are 
sufficiently small that the quintessence field (j)g is effectively 
frozen at the present epoch. When combined with axion dark 
matter and Affleck-Dine baryogenesis, this can provide an ex- 
plicit realization of the class of models we are considering 
here. 

The paper is organized as follows. In Section 2 we de- 
rive the general probability distribution for the ratio of the 
baryon density to the dark energy density at initial collapse, 
rg, and the baryon-to-dark matter ratio, r, for galaxies with 
equal baryon density at initial collapse. We also derive the 
constraints on possible models from the requirement that the 
probability distribution is well-defined. In Section 3 we ap- 
ply this distribution to a range of models and show that val- 
ues of rg and r which are consistent with the observed Uni- 
verse are quite natural in a number of models. In Section 
4 we present our conclusions. In the Appendix we discuss 
frozen quintessence models and show that they can consis- 
tently generate a domain-dependent frozen dark energy den- 
sity via quantum fluctuations. 



angles are equally probable 1 . However, we may consider 8 ( / 
and Qt, to be general random variables. The anthropic selec- 
tion condition that galaxies with equal baryon density at initial 
collapse are anthropically equally probable independently of 
r was then applied. Although an idealization, this condition is 
likely to be a good approximation for the range of baryon-to- 
dark matter ratio favoured by the resulting probability distri- 
bution, 0. 1 < r < 1 lfl2tl . It should be emphasized that the only 
quantities which vary between the domains are the densities. 
In particular, the primordial density perturbation is assumed 
to be the same in all domains. 

Under these condition, the probability density function 
f r {r) of the baryon-to-dark matter ratio, r, was shown to be 
determined purely by the statistical probability of the corre- 
sponding values of Q c / and 8^ 1 12], 
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where N is a normalization constant which depends on a and 
(3. Rather remarkably, this is completely independent of the 
specific value of the baryon density at initial collapse, p/, g . 

Physically, this distribution arises because different values 
of r can result in the same baryon density at initial collapse of 
galaxy perturbations. For example, if we increase the average 
dark matter density at a fixed temperature, the effect will be to 
increase the temperature at matter-radiation equality, resulting 
in earlier perturbation growth and so earlier initial collapse at 
a higher temperature. This increases the baryon density in 
galaxies at initial collapse. This can be compensated for by 
reducing the average baryon density at a fixed temperature, 
so keeping the baryon density at initial collapse fixed. The 
probability of r is then determined by the number of domains 
of different average dark matter and baryon density which can 
produce a given baryon density at initial collapse. 

We now consider the effect of a domain-dependent dark en- 
ergy density. We will consider a constant dark energy den- 
sity pg which depends on a random variable (j)g according 
to pg oc ty< where all values of ())g are considered equally 
probable. In the frozen quintessence models discussed in the 
Appendix, (|)g is the quintessence field and its potential is 

v(4> e ) oc 

The effect of increasing the dark energy density will be to 
damp perturbation growth once the total matter density is less 
than the dark energy density. Therefore initial collapse of the 
galaxy-forming perturbations must occur before the total mat- 
ter density Ptot is less than the dark energy density. In the fol- 
lowing we will consider a simplified model where the growth 



II. PROBABILITY DISTRIBUTION FOR r Q AND r 

In lfl2ll we considered dark matter and baryon densities 
which depend on random variables Q c i and 8/, according to 
p c i 8" and 8^. In the model of lfl2ll . 8^ and 8^ are the 
phases of massless fields which take random values during in- 
flation due to quantum fluctuations, such that all values of the 



We should make clear what we mean by a "domain". Quantum fluctuations 
randomize the fields at a number of e-foldings much larger than N ~ 60, 
when the observed Universe exits the horizon. At N ~ 60, there will then 
be approximately constant random values of the fields in any horizon-sized 
volume at N ~ 60, due to the random walk sum of fluctuations of wave- 
length much larger than the horizon at N t« 60 II ill . There will also be 
isocurvature fluctuations of wavelength close to the horizon atN~ 60, but 
these are small perturbations of the otherwise constant value of the random 
fields over the horizon. 
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of perturbations stops once pg > Ptot and is unaffected when 
Pg < Ptot- 

We next compute the probability for a domain to have a 
values of (rg,r) in the range (rg i,r\) to [tq^jj). Here rg = 
Pbg/p<2 is the ratio of the baryon density at initial collapse to 
the constant dark energy density. 

There are two cases. If rg > 1 then Ptot > Pg for all r 
at initial collapse, therefore the dark energy density has no 
effect on galaxy formation. The probability density / of the 
domains is then simply the product of the probability density 
f r of r and the probability density f rQ of rg, where the factors 
are independent 



integral becomes 
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The probability density for rg at a fixed value of pi, g is ob- 
tained from the flat distribution for (|)g (i.e. all (j)g having equal 
probability), 

Prob.oc^goc-^^ocJ-r- (3) 



On the other hand, if rg < 1 then, since p^ „ < pg, we require a 
sufficiently large dark matter density in order to have Ptot > 
Pg at initial collapse. This imposes an upper bound, r c , on r 
as a function of rg, 
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where r c = °° once rg > 1 . This in turn reduces the total num- 
ber of domains available at a given rg. 

The probability P to be in the range rg i to rg 2 and r\ to r2 
is then 

rr Q2 rMm( n .r c {r Q )) \ / r \ \ (o+h) 1 

"= '" ' ' 7¥{—r) 

r Q 

(5) 
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where the overall normalization N satisfies 
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P simply counts the total number of available domains with 
values of r and rg within the given ranges, subject to the con- 
ditions that pb g is constant and Ptot > Pg at initial collapse. 
An important feature is that the probability distribution for rg 
and r is independent of the specific value of the baryon den- 
sity at initial collapse pi g , which is assumed to be determined 
by anthropic selection. 

Eq. © is a well-defined probability as long as the integral 
in Eq. © is finite. As shown in 11211 . the integral over r is 
finite provided that a > P/3. The integral over rg can diverge 
at small rg. In the limit where the upper limit on the rg inte- 
gral is rg i <C 1, 'e(rg) is approximately equal to rg and the 
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Thus the condition for convergence as the upper limit rg i —> 



_3 1 (n + 1) 

4p 4a n 



1 >0 
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If the integral does not converge then the probability distribu- 
tion will strongly prefer small values of rg and r, correspond- 
ing to very large dark energy and dark matter densities relative 
to the baryon density. We expect the dark matter density will 
be anthropically cut-off when it becomes sufficiently large, 
because the increase in the dark matter density in the galaxy 
halo for a given baryon density will cause the radius of the 
disk to become smaller for a given disk angular momentum. 
In this case, assuming that the effect of increasing the dark 
matter density is dominated by the increase in mass of the 
halo (and not by an increase in the angular momentum of the 
disk), the density of stars in the disk will eventually become 
large enough to affect stable solar system formation. How- 
ever, we expect this to become important only at very small 
values of r and therefore not to alter the conclusion that very 
large dark matter and dark energy densities will be favoured in 
this case. Therefore it is unlikely that models with divergent 
probability distributions will be consistent with the observed 
Universe. This is a significant constraint on possible models. 



III. PROBABILITY OF THE DENSITY RATIOS IN 
SPECIFIC MODELS 



In this section we will consider some plausible examples 
of models, which are defined by a, P and n. We will con- 
sider models with n = 2 and n = 4. These can be realized 
by frozen quintessence models with and tyg potentials re- 
spectively (Appendix). For a and p we will consider models 
with (oc,p) = (1,1), (2,1), (1,2) and (2,2). (a,P) = (2,1) can 
be realized by axion dark matter with an approximately a 2 
potential combined with a baryon asymmetry determined by 
a domain-dependent CP-violating phase, as in Affleck-Dine 
baryogenesis 111 ill , (a, p) = (1,1) could be realized by asym- 
metric dark matter proportional to a domain-dependent CP- 
violating phase combined with a baryon asymmetry which is 
also proportional to a domain-dependent CP-violating phase. 

The range relevant to our domain is rg = 100 — 1000 and 
r = 0. 1 — 1 . The range of r is determined by the magnitude of 
the observed baryon-to-dark matter ratio, r = 0.2. To estimate 
the range of rg which is consistent with the observed Uni- 
verse, we need to estimate the dark energy density when the 
galaxy perturbations first go non-linear and begin to collapse. 
The red-shift at the initial collapse of a galaxy perturbation, 
Zcolli is related to rg and the present density of baryons and 
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dark energy by 



rg = (1 +z C oll)' 



(9) 



With £l b = 0.04 and Q. Q = 0.73 we find that rg w 95 if z co ll ~ 
11 and rg « 1000 if z co « « 25. Therefore rg in our domain 
is in the range 100-1000 if z co u is in the range 11 — 25. This 
range is consistent with estimates of the time of formation of 
proto-galaxies. 

We first consider which models have convergent probability 
distributions. For n = 2 the convergence condition is 



3 1 „ 
p-« >2 - 
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This can only be satisfied by (a, p) = (2, 1). For n = 4 the 
convergence condition is 



3 1 , 



(ii) 



This can be satisfied by (a,p) = (2,1) and (1, 1). Models 
with (a, p) = (1,2) and (2,2) are generally ruled out by con- 
vergence when n = 2 or 4. 

We next discuss the probability distributions for the three 
convergent models and their compatibility with the observed 
Universe. 

We first consider the case n = 2 and (a, p) = (2,1). In Ta- 
ble 1 we give the probability for rg and r to be within different 
ranges of values, where we express the probability in terms of 
the percentage of the total number of domains. This shows 
that the probability of rg and r being within the observed 
range in this model is small, with only 0.9% of domains hav- 
ing rg = 100 — 1000 and r = 0.1 — 1. This can be compared 
with the most likely range, rg = 0. 1 — 1 and r = 0.1 — 1, which 
has 13.5% of the domains and is therefore 14.8 times more 
likely. Thus domains with dark energy densities which are 
much larger than those we observe are strongly preferred in 
this model. This is also apparent from Figure 1, which shows 
the strong preference for domains with much smaller rg (i.e. 
larger dark energy density relative to baryons) than we ob- 
serve. The prospects for the observed value of r are better. 
Figure 2 shows that values of r in the range 0.1-1 are gener- 
ally the most likely or close to the most likely when rg > 0. 1 . 

Consistency with the observed Universe is much improved 
in the model with n = 4 and (a, p) = (2, 1). From Table 2 we 
see that the probabilities for all ranges of rg and r are not very 
much different. The probability (12.1%) of the most likely 
range of rg and r (rg = 1 — 10 and r = 1 — 10), is only 3.5 
times larger than the probability (3.4%) of range correspond- 
ing to our domain (rg = 100 — 1000 and r = 0.1 — 1). This is 
also clear from Figures 3 and 4. From Figure 3 we see that the 
probability distribution of rg is more evenly spread across the 
ranges of rg than in the case n = 2, with the peak of the distri- 
bution being at larger values of rg. Therefore values of rg in 
the range consistent with the observed Universe are a natural 
possibility in this model, with a probability which is not very 
much smaller than the other ranges of rg. From Figure 4 we 



see that the probability of r in the range 0. 1 - 1 .0 is large for the 
most likely ranges of rg. Therefore the model with n = 4 and 
(a, p) = (2,1) can naturally account for the observed coinci- 
dence of the dark energy, dark matter and baryon densities. 

The case n = 4 and (a, p) = (1,1) is similar to n = 4 and 
(a, p) = (2,1). From Table 3 and Figures 5 and 6 we see that 
the distribution of values of rg is similar to but slightly less 
strongly peaked than the case (a, P) = (2, 1). The probability 
of rg and r being in the range consistent with our observed 
Universe is again relatively large (3.6%) compared with the 
most likely range (rg = 0.1 — 1 and r = 0.1 — 1), which is 3.1 
times more likely with 11.1% of the domains. Therefore it 
is again quite reasonable in this model to find ourselves in a 
domain with the observed dark energy, dark matter and baryon 
densities. 

The preference of the observed large value of rg in our do- 
main for models with larger values of n can be easily under- 
stood. The gain in probability in going to larger (j)g is smaller 
in this case, since an increase in (j)g produces a larger increase 
in the dark energy density pg °< (j)g. Therefore larger dark en- 
ergy densities (corresponding to smaller rg) are less strongly 
favoured. 
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1-10 


10-100 


0.001-0.01 


9.8 














0.01-0.1 


8.0 


9.6 


0.03 








0.1-1 


2.5 


9.7 


13.5 


2.8 


0.15 


1-10 


0.8 


3.1 


8.9 


10.2 


4.6 


10-100 


0.26 


1.0 


2.9 


3.2 


1.4 


100-1000 


0.08 


0.31 


0.91 


1.0 


0.46 


1000-10000 


0.04 


0.11 


0.33 


0.37 


0.17 



TABLE I: Percentage of domains within ranges of (rg,r) for n = 2 and 
(oc,P) = (2,l). 
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0.01-0.1 


3.0 


4.2 


0.02 








0.1-1 


1.7 


6.4 


10.0 


4.0 


0.14 


1-10 


0.95 


3.6 


10.6 


12.1 


5.3 


10-100 


0.54 


2.1 


6.1 


6.9 


3.0 


100-1000 


0.30 


1.2 


3.4 


3.9 


1.7 


1000-10000 


0.20 


0.77 


2.3 


2.6 


1.2 



TABLE II: Percentage of domains within ranges of (rg,r) for n = 4 and 
(o,P) = (2,l). 



IV. CONCLUSIONS 

We have shown that is possible to explain the coincidence 
of the observed dark energy, dark matter and baryon densities 
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FIG. 1: Percentage of domains as a function of the baryon-to-dark 
energy ratio at initial collapse tq for n = 2 and (a, (3) = (2, 1) 




FIG. 2: Percentage of domains as a function of the baryon-to-dark 
energy ratio at initial collapse rg and the baryon-to-dark matter ratio 
rforn = 2and(a,P) = (2,l) 
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FIG. 3: Percentage of domains as a function of the baryon-to-dark 
energy ratio at initial collapse yq for n = 4 and (a, (3) = (2, 1) 
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9.9 


11.1 
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11.1 
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1.1 


10-100 


1.1 


3.2 


6.3 


3.8 


0.65 


100-1000 


0.65 


1.8 


3.6 


2.2 


0.36 


1000-10000 


0.43 


1.2 


2.4 


1.4 


0.24 



TABLE III: Percentage of domains within ranges of (rg,r) for n = 4 and 
(a,p) = (l,l). 




FIG. 4: Percentage of domains as a function of the baryon-to-dark 
energy ratio at initial collapse rQ and the baryon-to-dark matter ratio 
rforn = 4and (a, (3) = (2,1) 



via domains of varying densities generated by random vari- 
ables combined with anthropic selection which is primarily 
sensitive to the baryon density at the initial collapse of galaxy- 
forming perturbations. 

We have derived the complete probability distribution for 
the number of domains with fixed baryon density at initial 
collapse as a function of the baryon-to-dark energy ratio at 
initial collapse, rg, the baryon-to-dark matter ratio, r, and the 
model parameters n, a and p, where the respective densities 
depend on the random variables (j)g, Q ( / and % according to 

PQ x §Q' Pdm x 9™ and p& « 0^. This allows us to com- 
pute the probability of finding ourselves in a domain with a 
given magnitude of rg and r. The probability distribution for 
rg and r is essentially statistical in nature, being independent 
of the anthropically-selected baryon density. A natural inter- 
pretation of the random variables is that they are effectively 
massless fields during inflation, whose domain-dependence is 
generated by quantum fluctuations. 

An important constraint on possible models is the require- 
ment that the total probability converges to a finite value. We 
find that for n = 2, which could be realized by a frozen 
quintessence model, this is only possible for models with 
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FIG. 5: Percentage of domains as a function of the bary on-to-dark 
energy ratio at initial collapse yq for n = 4 and (a, (3) = (1,1) 



r 




FIG. 6: Percentage of domains as a function of the baryon-to-dark 
energy ratio at initial collapse rg and the baryon-to-dark matter ratio 
rforn = 4 and (oc,P) = (1,1) 



(a, p) = (2, 1), while for n = 4, which could be realized by 
a frozen quintessence model, it is possible for models with 
(a,P) = (2,1) and (1,1). (oc,p) = (2,1) could be realized 
by, for example, axion dark matter combined with Affleck- 
Dine baryogenesis (or, more generally, any form of baryo- 
genesis based on a domain-dependent CP-violating phase), 
while (oc,p) = (1,1) could be realized by asymmetries ofboth 
dark matter and baryons which are determined by domain- 
dependent CP- violating phases. 

Comparing the probability for rg and r with the ranges con- 
sistent with the observed Universe, we find that the model 
with n = 2 is not favoured, with only a small fraction of the 
domains being consistent with the observed Universe, making 
our domain atypical. On the other hand, we find that for mod- 
els with « = 4, our domain is a more natural possibility. The 
distribution of dark energy densities has broadly similar prob- 
abilities over a wide range of rg, while for larger values of rg 



the most likely value of r lies in the range 0.1-10. The do- 
mains with rg and r in the range consistent with the observed 
Universe are fairly typical, with only about 3.5 times smaller 
probability than the most likely range of values of rg and r. 
Thus it is quite natural to find values of the dark matter and 
dark energy densities relative to the baryon density which are 
consistent with those we observe, even though larger values 
of the dark energy density are more probable. 

The fundamental assumption of the model is that anthropic 
selection is primarily sensitive to the baryon density at ini- 
tial collapse and is relatively insensitive to the value of r. 
For baryon densities at initial collapse similar to the observed 
Universe, anthropic sensitivity to r is likely to become impor- 
tant only at very large dark matter densities, corresponding to 
small values of r 0.1. These values of r make only a small 
contribution to the probability distribution for r, therefore the 
model should be reliable in this case. Thus it should be pos- 
sible to explain the dark matter and dark energy densities in 
our domain via domain-dependent densities. Nevertheless, it 
would be interesting to consider the effect on the final state 
of galaxies of varying the baryon density and baryon-to-dark 
matter ratio at initial collapse. This would allow us to assess 
the anthropic selection of the baryon density and the likeli- 
hood of the observed baryon density, and to quantify the an- 
thropic sensitivity to r. This is quite challenging, as we would 
need to understand the effect of different initial baryon and 
dark matter densities on the typical size of galaxies, the pro- 
cess of gas cloud collapse and star formation and the final state 
of the galactic disks. 

In conclusion, we find that it is possible to naturally account 
for the coincidence of dark energy, dark matter and baryons in 
a framework in which all densities are determined by random 
variables and the dominant anthropic selection parameter is 
the baryon density in galaxies at the initial collapse of galaxy- 
forming perturbations. 
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Appendix: Frozen quintessence as a model for a 
domain-dependent dark energy density 

As an example of a model of a domain-dependent dark en- 
ergy density, we will consider a simple frozen quintessence 
model based on a potential V((j)g) = K§q , where the effective 

mass squared V "((j)g) of (j)g is much less that the present value 
of the Hubble parameter squared, H} } . (Although V "((j)g) 
could be close to Hg, this would amount to an unlikely co- 
incidence, therefore we will consider V ((|)g) <C H% to be the 
natural condition.) The domains will then be generated by 
quantum fluctuations of the effectively massless quintessence 
field. As a background inflation model we will consider (j) 2 
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chaotic inflation. 

The condition V (<|)g) <C H% at present imposes an upper 
bound on the value of K, 



K < K m 



(A-l) 



Q 



The condition for the energy density to be dominated by dark 
energy at present is that H 2 = V((j)g) /3M 2 . (M p is the reduced 
Planck mass throughout, 2.4 x 10 18 GeV.) Therefore 



L/« 



(A-2) 



Combining Eq. ( lA-U and Eq. ( 1A-2I ) implies that 



^) e »(3«(n-l)) 1 / 2 M p . (A-3) 

Therefore (|)g ^ M p is necessary to have frozen quintessence. 

We next show that must be generated by quantum fluc- 
tuations in order to satisfy (])g 3> M p . We assume that is 
effectively massless during inflation. We also assume that the 
length scale of initial fluctuations of the fields is characterised 
by the initial chaotic inflation horizon. In this case the initial 
fluctuations 8(j)g cannot be larger than M p . This follows from 
the contribution of the gradient term to the energy density 



Pe ~ ( Vc M 2 



n 2 m 



(A-4) 



Therefore pg < p = 3M 2 H 2 implies that 8(j)g < M p . Therefore 
if we wish to have (|)g 3> M p in domains atN — 60, we require 
that quantum fluctuations of (j)g can subsequently drive the 
field to large values. 

We next show that this can easily be achieved in (j) 2 chaotic 
inflation. We consider the initial condition to be H ps M p , 
corresponding to a Planck-scale initial energy density, as ex- 
pected from chaotic initial conditions. The initial energy den- 
sity of the inflating patch is V(<j>/), where (]), is the initial value 
of the inflaton field. Therefore 



V6Mj 

fo« 

mis/ 



(A-5) 



The number of e-foldings until the end of inflation N for a 
given (]) is 



N = 2kG^ 2 -^,J^2kG^ , 



(A-6) 



where § em ( ps \f2M p (<^ (]),■) from the condition T\(^ en j) ps 1. 
Therefore there is a very large total number of e-foldings from 

(j) = (j); to § end , 

, 3 fM„\ 2 
2 \m<j, / 

»u ps 1 .5 x 10 13 GeV is necessary to generate the correct mag- 
nitude of primordial density perturbation, therefore Ntot ~ 
4x 10 10 . 



The average quantum fluctuation of the quintessence field 
on the scale of the horizon per e-folding is 8(|)g = H/2n. If 
H is constant, the resulting random walk of the field will re- 
sult in a Gaussian distribution for the field in a given horizon 
volume. In fact, H is not constant during (|) 2 chaotic inflation, 
since H — (2m^/3) 1 / 2 at N e-foldings. Therefore the proba- 
bility distribution for the effectively constant value of (j)g in a 
domain at N ps 60 will not be exactly Gaussian, since the ran- 
dom walk step of the (j)g field is time-dependent. Nevertheless, 
we can compute the variance a of the field by summing the 
variance due to small numbers of e-foldings, since H varies 
on a time scale much smaller than an e-folding. Therefore 



a 2 =< m >= y o 



Ntot H 2 (N) H Z (N- 
dN = 



4n 2 



4k 2 



tot) Ntot 
- x 



(A-8) 

With Ntot ~ 4 x 10 10 we obtain o w 2 x 10 4 M p . 

Thus quantum fluctuations of the quintessence field can 
generate superhorizon-sized domains at = 60 with effec- 
tively constant values of (j)g which are much larger than M p . 
For values (|)g <C O, all values of (j)g will be equally probable, 
as in the case of an exact Gaussian distribution. This can be 
seen since most of the variance is generated at N > Ntot/2, 
where H is approximately constant and the distribution of (j)g 
is close to Gaussian. 

(j)g at present can remain constant at all values up to the 
variance a if v"(g) -C H 2 . This is true if 



n(n - l)a"~ 



(A-9) 



For n = 2 this requires that mg <C H , since K = /2. There- 
fore niQ <C 10~ 42 GeV is necessary. For n = 4 the condition 
becomes 



K< 



2k 2 H 2 
3N T otM 2 



10 



-130 



(A-10) 



Therefore a severe tuning of the self-coupling is necessary 
to allow (j)g to take any value out to the variance. The value of 
())g in our domain could be much smaller than the variance, in 
which case larger K could be considered. However, even with 
the smallest possible value of (|)g, (j)g ps M p , the n = 4 model 
still requires K < 10~ 120 . We note that a sufficiently large 
maximum value of (|)g compared with M p is necessary to al- 
low a large enough range of dark energy density for anthropic 
selection to function. 

A combination of super-Planckian (j)g and extremely small 
mass and self-coupling terms for (j)g is therefore necessary to 
have a successful frozen quintessence model with dark en- 
ergy domains generated by quantum fluctuations. This ap- 
pears to be a rather extreme set of assumptions. However, 
the assumptions of the frozen quintessence model are sim- 
ply a more extreme version of the assumptions behind the 
(j) 2 chaotic inflation model. In the (|) 2 chaotic inflation model 
§(N) ss 2VNM p > M p at N e-foldings and m$/M p ps 10~ 5 . 
Therefore the difference between the tuning of the two mod- 
els is only a question of degree, which will depend on the 
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measure for the parameters (the difference in tuning appears 
less extreme if viewed logarithmically). 

Finally, it is important to check that quantum fluctuations 
of generated at N w 60 do not lead to dangerous fluctua- 
tions of the dark energy on observable scales. If we consider a 
fluctuation 8(])g then the energy density in the fluctuation once 
within the horizon is 



Thus 



Sp e « V ((j) e )8(j) e , 



(A-ll) 



where we assume that the energy of the fluctuation is dom- 
inated by its potential, which is true if k 2 8(j)g < V ((j)g)8(j)g, 
where k is the wavenumber of the fluctuation. If we consider 
a fluctuation entering the horizon at the present time, corre- 
sponding to a perturbation which exited the horizon at N « 60, 
then 8(j) e w H(N = 60) /2jc. With V(<S?q) « §" Q , we find 



8p e _ n8$Q _ nH(N = 60) 



(A-12) 



PQ 



(A-13) 



Since typically (j)g > M p , it follows that Spg/pg <C 10~ 5 
and so quantum fluctuations of the quintessence field will 
have a completely negligible effect. The condition k 2 8(j)g < 

V'($q)8$q is satisfied in the case k 2 ~ Hg if H(N = 60) < 
M 2 /(j)g, which is generally true for the models considered 
here. 

Thus frozen quintessence with a (|)g potential can consis- 
tently generate domains of different dark energy density. The 
probability distribution of is flat, so (j)g acts as a random 
variable which generates an effective cosmological constant 
in each superhorizon-sized domain at present. This provides 
an existence proof for such models, although, as with all 
quintessence models, rather extreme assumptions the regard- 
ing mass and coupling terms are necessary. 



[1] S. Weinberg, Phys. Rev. Lett. 59, 2607 (1987). 

[2] S. Weinberg, Rev. Mod. Phys. 61 (1989)_L 

[3] J. McDonald, JCAP 0701 (2007) 001 |hep-ph/0609T26). 

[4] J. McDonald, Phys. Rev. D 83, 083509 (2011) 

larXiv: 1009.32271 [hep-ph] ] . 
[5] J. McDonald, Phys. Rev. D 84, 103514 (2011) 

larXiv: 1108.46531 [hep-ph] ] . 
[6] Y. Cui, L. Randa ll and B. Shuve, JHEP 1204, 075 (2012) 

larXiv: 11 12.27041 [hep-ph]]. 

[7] Y. Cui and R. Sundrum, arXi v:1212.2973l [hep-phi. 



[8] S. Davidson and M. Elmer, JHEP 1210, 148 (2012) 

larXiv: 1208.055 II [hep-ph] ] . 
[9] A. D. Linde, Phys. Lett. B 201, 437 (1988). 

[10] A. D. Linde, Phys. Lett. B 160 (1985) 243^ 

[11] J. McDonald, JCAP 1210, 005 (2012) larXiv:1207.2135i rhep- 

ph]]. 

[12] J. McDonald, larXiv: 1211 .44721 rhep-phl . 
[13] R. R. Caldw ell and E. V. Linder, Phys. Rev. Lett. 95 (2005) 
141301 |astro-ph/0505494) . 



